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Abstract 

Within a microscopic approach which takes into account RPA configurations, the 
single-particle continuum and more complex \plh®phonon configurations isoscalar 
and isovector Ml excitations for the unstable nuclei 56 > 78 Ni and 100 > 132 Sn are cal- 
culated. For comparison, the experimentally known Ml excitations in 40 Ca and 
208 Pb have also been calculated. In the latter nuclei good agreement in the centroid 
energy, the total transition strength and the resonance width is obtained. With 
the same parameters we predict the magnetic excitations for the unstable nuclei. 
The strength is sufficiently concentrated to be measurable in radioactive beam ex- 
periments. New features are found for the very neutron rich nucleus 78 Ni and the 
neutron deficient nucleus 100 Sn. 
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1 Introduction 



There is a rich history in the study of Ml excitations in nuclei (see ref. || for a recent 
review). For more than two decades the experimental effort has concentrated on stable 
magic nuclei . Unfortunately, in the valley of stability, the number of cases is quite 
limited which rules out systematics of the quality we are used to from electric resonances 
||. With the availability of radioactive beam facilities this may change since far from 
stability new nuclei become available with similar shell characteristics as the stable magic 
ones ||. Aside from systematics such nuclei are also of great astrophysical interest (see 
references in 

The aim of the present paper is to extend the microscopic models which describe 
well the known resonances to the new region accessible by radioactive beams. As is well 
known, the RPA with phenomenological Landau-Migdal interactions is able to reproduce 
the excitation energies of known isoscalar and isovector Ml excitations with a universal 
set of spin-interaction parameters g and g' . Including effective spin g-factors ~ 0.7 — 0.8 
the observed total transition strengths are also reproduced || |5], [7|, BJ. The RPA fails, 
however, to account for the resonance width which, especially in heavy nuclei, is caused by 
a coupling to more complex configurations. For the Ml case there are several approaches 
for including such configurations || 10, ll|. Concerning the quenching of Ml strength the 
most reliable approach seems to be the use of phenomenological values for the g-factors. 



By now, these are very well establish from a variety of experiments || [12| and are universal 
for all nuclei measured. For predictions in unstable nuclei knowledge of global interaction 
parameters and g-factors is of crucial importance. 

Most realistically, the more complex configurations are described as one-particle one- 
hole states coupled to low-lying surface vibrations {Iplh <g> phonon states) using known 
parameters. Recently a formalism for including such configurations in conjunction with 
the single-particle continuum has been successfully used in the description of giant electric 
resonances ]T3|, [14[]. In particular, the widths are reproduced satisfactorily. We shall 
employ this approach here in the study of Ml transitions in unstable nuclei. 

The theory will be briefly described in sect. 2. In sections 3. and 4. calculations 
for the known Ml states in 48 Ca and 208 Pb are presented which show the quality of the 
theoretical results. The same interaction parameters, effective charges and a standard 
mean field of the Woods-Saxon type are then used to predict spectra for unstable nuclei. 
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2 Theory 



The approach which we shall use has been developed in the framework of the consis- 
tent Green's function method and is based on Migdal's Theory of Finite Fermi Systems 
(TFFS) [jn|. The main physical idea is to include lplh ® phonon configurations instead 
of "pure" 2p2h ones and thus to make use of the fact that, in magic nuclei, the squared 
phonon creation amplitude is a small parameter g 2 < 1. This enables one to restrict 
the calculation to particle-phonon coupling terms of order g 2 . Furthermore only a small 
number of phonons of maximum amplitude (maximal g 2 ) need to be considered, i.e. the 
most collective low-lying ones. This makes it possible to use known parameters of the 
TFFS which determine the local effective interaction and quasiparticle charges (taking 
phonons into account explicitly gives non-local contributions). In addition, a small num- 
ber of phonons greatly reduces the computational effort. This physical assumption has 



been confirmed in several calculations of Ml excitations [11 , the giant dipole resonance 
m |TJ as well as E2, EO resonances |T§ . 



The theory is formulated, most efficiently, in coordinate space. For an applied external 
field V° with frequency oj, the change in density 5p = p — p from equilibrium is given by 

6p(r,u) = -J dv'A{v,v',u){e q V Q {v')-T{v')5p{v',u)) (1) 

and the frequency distribution determines the excitation spectrum of the system. Here e q 
denotes the local quasiparticle charge and T the quasiparticle interaction. As discussed 
above, these can be taken from the TFFS. The generalized propagator A contains the 
RPA part as well as the lplh <E> phonon configurations and was derived in ref. |17| in a 
representation of the discrete single-particle basis {<fi\,£\}. The tilde indicates that these 
states have been corrected for mean field contributions from the particle-phonon coupling, 
already included in phenomenological potentials of the Woods-Saxon type. Since these 
contributions are included explicitly they have to be removed from the single-particle 
potential to avoid double counting. To include single-particle emission, the discrete prop- 
agator v4 12 34 has to be augmented by a continuum part and takes the form 

A(r,r» = A™{tJ,lj) + £[A 1234 H - if^M^^^W^W^^^Kr') (2) 

1234 



where A^f is the continuum RPA propagator (for details see ref. |18|, [19] 
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The Ml transition probabilities between the ground- and excited states are determined 
by the strength function 

^^ = EK-Ie^°|0)| 2 ^-^). (3) 
Via the optical theorem this is related to Sp as 

dB [ M ^ =-hm( dre q V°(r)6p(T, u) (4) 

and the solutions of Eq. ([!]) determine dB(Ml) / duj. 

The summation over single-particle states in Eq. (0) is usually performed for two shell 
above and below the Fermi level. Aside from providing a finite width above the continuum 
threshold our method has the distinct advantage, over discretized calculations of this type, 
that the matrix dimension still remains manageable. 

A major difference between the approach presented above and similar approaches 
20| , |21~1 , |22|| is the consistent inclusion of ground state correlations (GSC) beyond the RPA. 



There are two kinds of such 2p2h (to be more exact lplh®phonon in our case) correlations. 
One is GSCs without "backward going" diagrams containing the quasiparticle-phonon 
interaction. The second one which includes these diagrams has new poles which generate 
new states. In the RPA, GSCs are caused by non-pole diagrams, only, which cannot 
produce additional excited states. In this sense our theoretical extension is qualitatively 
different from the RPA. For Ml excitations in magic nuclei the consequences have been 
discussed in Refs. |TT|, 1231 without taking into account the single-particle continuum. For 



the spin-saturated nuclei 16 and 40 Ca, where strong Ml excitations are not allowed 
because of missing spin-orbit partners, one obtains Ml excitations which describe the 
experimental data reasonably well. 

3 Approximations and Parameters 



The main physical approximations have been described above and in [jl3| , [Ljfl . Here we 
present the parameters used for the present calculations. 

According to the selection rules only the spin-dependent part of the Landau-Migdal 
interaction 

T = C (g + g'r ■ r')a ■ o^n - r 2 ) (5) 
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enters in the Ml excitations. In order to predict the Ml states in the unstable nuclei 
reliably we have slightly adjusted the g' parameter to obtain the well-known 1 + levels in 
48 Ca. This gives g' = 0.86 instead of g' = 0.96 used previously § [II], |T7[ (C„ = 300 
MeV fm 3 ). For the parameter g we have used g = —0.05. 
The local quasiparticle charge was determined from [15 



epy* = (l-^ + Kl-Uf + ^ + l^-l/^ 

e n V on = + [(1 - e.) 7 » + e.7 p - (6) 
where 7 P = 2.79/zo, 7™ = — 1.91/xo, A*q = eh/2m p c and 

e = c = o.i, er = = -0.03 (7) 

as obtained earlier [8], [IT], |23], I7|| . These values yield for the spin local charge e p q = 



0.647 p and e™ = 0.747™. 

The low-lying collective phonons (listed in Table 1 for the various nuclei considered) 
which are taken into account in the lplh <S> phonon coupling have been calculated within 
RPA using the following known Landau-Migdal parameters 

f m = -0.002, f' ex = 2.30, fi = 0.76 
g = -0.05, g' = 0.96, C = 300MeVfm 3 (8) 

for all nuclei. Since the spin-dependent contributions to the low-lying phonon spectrum is 
small, the difference between g' = 0.86 and g' = 0.96 has no noticeable consequences for 
the results. The parameters f ex has been fitted to available phonon energies and ranges 
from -5.0 to -3.5 which is not far from the value f ex = —3.74 used earlier [ 151 , |T4| , 



The single-particle states have been calculated from a standard Woods-Saxon potential 



2q| . To get good agreement with measured single-particle energies which are available for 
48 Ca, 208 Pb, 56 Ni |g, 132 Sn ^ || and, to some extent, for 100 Sn [||, || the well depth 



of the central part of the potential has been adjusted, changing the depth parameter by 
less than 5%. A list of single-particle energies, thus obtained, is given in Table 2. Another 
check is the energy of the gg/2 (p) — > ^(n) Gamow- Teller transition in 100 Sn. We find a 



transition energy of 7.6 MeV which coincides with the results of Ref. [p8 |. 

As mentioned above, the coupling of single-particle states to phonons gives a contribu- 
tion to the observed single-particle energies. When adjusting the Woods-Saxon potential 
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to empirical energies such contributions are already accounted for. In order to avoid 
double counting the particle-phonon contributions should be removed. This leads to a 
"refined" basis {<fi\, £\}. The corrected energies i\ are obtained by solving the non-linear 
equations 

e A = ex- M\\(t\) 

■ = fr-^EKAWi'i , 1 ~- n " , + -* , } o» 

2 3 + 1 t£' le x -e\i -uj s e x - e x < + uj s J 

where A = {n, 

Finally, in order to account for instrumental resolution in the experiment we have 
included an energy averaging parameter A. 

4 Results and Discussion 

The results of our calculations are presented in Table 3 and in Figs. 1-10. The continuum 
RPA results are denoted by lplh+ continuum while lplh+2p2h+continuum indicates the 
results including the lplh®phonon contributions. To elucidate the role of the additional 
ground state correlations we quote results with and without (gs(±)) in the figures. The 
mean energies listed in Table 3 are defined as 

= Ei Bi(Ml T) (10) 

Fig. 1 shows the calculated strength distribution for 48 Ca. In the continuum RPA (dotted 
line) this is dominated by the If 7/2(71) — > 1/5/5(71) transition. After inclusion of more 
complicated configurations the measured centroid energy of 10.23 MeV is quite well re- 
produced (E th = 10.36 MeV) while the calculated transition strength is somewhat larger 
than in the experiment. The additional ground state correlations have no noticeable effect. 
In 208 Pb the two particle-hole transitions 1/111/2(2?) — * 1^9/2(2?) arid l«i3/2(n.) — > lin/2(n) 
form an isoscalar and isovector transition within the RPA. Taking into account the more 
complex configurations we obtain reasonable agreement with experiment for both the 
isoscalar 1 + level at E exp =5.85 MeV and for the isovector Ml resonance (Fig. 2). A 
rather strong quenching has been obtained so that our value J2B(M1 f) = 11.6/Xq is less 
than the experimental value of 15.6/Xq p9| - The calculated centroid energy (E th = 7.74 
MeV) is slightly above the measured one (E ex = 7.3 MeV). The width which can be 
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deduced from Fig. 2 is in good agreement with the experimental value of ~ 1 MeV. The 
effect of additional GSC is more pronounced than in 48 Ca: it increases the strength in the 
energy interval 9.9-15.5 MeV by about 3/ig (see Table 3). 

Table 3 and Fig. 3 give the results for 56 Ni M. As with 208 Pb this is a non-spin 
saturated nucleus with the I/7/2 shell filled for protons and neutrons. Thus an isoscalar 
and isovector state can be formed. The influence of the lplh <g> phonon states is quite 
noticeable reducing the strength in the peak by more than a factor of two. As can be 
seen from Fig. 3 and, in particular, from Fig. 8 (here the averaging parameter is 20 keV 
only) there is no fragmentation width for the resonance peak. The latter result also seems 
to hold for 132 Sn (Figs. 6 and 10). For 56 Ni we confirm the result of ref. [ 30 1 where no 
fragmentation width was found in a continuum RPA calculation. 

We now come to the "very exotic" nuclei 78 Ni and 100 Sn. In the former the 1/7/2(1?) 
and lgg/2 (n) are occupied while in the latter both 1*79/2 shells are full. As seen in Fig. 4 
78 Ni has an asymmetric resonance shape with a width of ~ 1 MeV. The Ml resonance 
in the very neutron deficient 100 Sn is split into two major peaks with E\ = 9.8 MeV and 
E2 = 10.5 MeV (Fig. 5). The low-lying isoscalar level is spit also. These features are 
caused by the inclusion of the lplh® phonon configurations and for 100 Sn also by the very 
small value of the proton binding energy B p = 2.91 MeV. 

5 Conclusion 

Within a microscopic approach which includes the continuum RPA as well as lplh <g) 
phonon configurations we have calculated isoscalar and isovector Ml excitations in the 
unstable nuclei 56,78 Ni and 100 ' 138 Sn. To judge the quality of our predictions we have 
also given results for the "known" Ml excitations in 48 Ca and 208 Pb. For the latter good 
agreement with experiment in the centroid energy, the transition strength and the width is 
obtained. This, together with the fact that we have used known universal parameters for 
the effective interaction and the local Ml charge, gives us confidence that the predictions 
for the unstable nuclei are realistic. 

The location of the isovector Ml resonances agrees rather well with the RPA calcula- 
tions. In all cases studied, we find a considerable influence of the particle-phonon coupling 
on the width and the strength distributions. As compared with the RPA there is notice- 
able strength in the high-energy tails and a decrease of strength in the resonance region. 
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The damping effects are particularly pronounced in the extremely neutron rich nucleus 
78 Ni as well as the neutron deficient 100 Sn. We did not find a simple A-dependence of the 
resonance widths. All these features depend delicately on the interplay between the shell 
structure and the low-lying vibrational modes. 

For approaches, like the one used here, which are based on a phenomenological single- 
particle scheme and the effective Landau-Migdal interaction it is important to experi- 
mentally confirm the results presented here. In particular, it would give reassurance of 
the universality (A-independence) of the Migdal parameters (||) and (0) and hence would 
allow to generalize the present approach to non-magic nuclei and to use it for many nuclei 
far from stability. 
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Table 1 : Characteristics of the low-lying phonons used in the calculations 



J" 


E, 
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e 2 fm 2L 


r 


E, 

Mev 


B(EL)T, 
e 2 fm 2L 


r 


E, 

MeV 
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e 2 fm 2i 


ioo Sn 






48 Ca 


2+ 


2.86 
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2+ 


4.06 


942 


2+ 


3.83 


81.6 


3" 


3.56 


5.33 10 3 


3" 


4.34 


5.0 10 4 


3" 


4.50 


1.12 10 4 


3- 


5.03 


5.99 10 3 


5" 


4.91 


3.46 10 6 
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Pb 


5- 


3.57 
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5- 


5.61 
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3" 
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4+ 
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2.33 10 3 
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1.68 10 5 
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6+ 
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4+ 


4.24 


1.80 10 6 


4+ 
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6+ 
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2.65 10 9 


6+ 


4.75 


1.14 10 9 


6+ 
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6+ 
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6+ 
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3.20 
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3.70 
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4+ 
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6+ 


4.00 
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6+ 


5.18 


1.36 10 8 


6+ 


4.94 


4.17 10 7 














3" 


4.38 


8.32 10 3 














5" 


5.04 


1.18 10 6 
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Table 2 : The single-particle levels for the unstable nuclei used in the 

calculations 
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Table 2 : The single-particle levels for the unstable nuclei used in the 

calculations 
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-9.87 
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2.95 
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3.85 
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2d3/2 
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5.85 
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3p3/2 
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lil3/2 






7.85 
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Table 3 : Characteristics of Ml excitations 
(E is ,E iv ,E are given in MeV,B(Ml) | and EB(M1) | in fi 2 ) 





lplft + cont. 


lplh + 2p2h + cont. 
without 2p2h GSC 




E- 


B(M1) T 

v. y i 


E- 


}^B(M1) T 
(interval) 




B(M1) T 


E 


(interval) 


48 Ca 






10.68 


8.64 






E- — 
10.35 


6.55 


Figs. 1,7 








(~ 10.68) 








(~ 10.35) 
8.24 
(9-19) 




5.63 


0.49 


7.97 


18.16 


5.72 


0.84 


7.66 


11.87 


208 Pfe 
Fig. 2 








(~ 7.97) 
(5-15.5) 








(6.3-8.7) 
(8.7-9.9) 

4.98 
(9.9-15.5) 

19.09 
(5.0-15.5) 




6.56 


0.31 


10.12 


11.31 


6.7 


0.24 


9.93 


4.8 


56 Ni 








(~ 10.12) 
11.43 








(8.5-11.0) 
5.58 


Figs. 3,8 








(5.5-15) 
11.73 
(6-45) 








(11-15) 
10.62 
(5.5-15) 




5.88 


0.66 


10.16 


16.0 


6.12 


0.64 


9.3 


10.2 


78 JVi 








(~ 10.16) 
16.7 








(8.5-10.2) 
6.8 


Figs. 4,9 








(6.5-15) 








(10.2-15.0) 
17.64 
(5.5-15) 




6.8 


1.5 


10.15 


14.1 


6.5 


0.3 


8.9 


1.3 


lm) Sn 








(~ 10.15) 








(7-9.4) 










15.6 


6.78 


1.4 


9.77 


6.6 


Fig. 5 








(6-13) 


(E = 
6.67) 




10.5 
(E = 
10.0) 


(9.4-10.15) 
3.2 

(10.15-10.7) 

3 

(10.7-13) 
15.8 
(6-13) 




5.8 


0.4 


8.78 


19.0 


5.8 


0.4 


8.45 


13.2 


132 Sn 








(~ 8.78) 
19.6 








(~ 8.45) 
6.0 


Figs. 6,10 








(5.5-15) 








(9-15) 
19.2 
(5.5-15) 



Table 3 : (continuation) 





Iplh + 2p2ft + cont. 
with 2p2h GSC 


Experiment 




Eis 


B(Afl) T 


E 


^B(Ml) T 
(interval) 




B(M1) T 


E iv or E 


B(M1) T or ^JB(MI) T 
(interval) 


48 Ca 
Figs. 1,7 






10.36 


6.12 
(~ 10.36) 
9.68 
(9-19) 






E — 
10.32 


3.9±0.3[31] 

5.3±0.6[31] 
(7.7-12.7) 


208p fe 

Fig. 2 


5.72 


0.84 


7.74 


11.57 
(6.3-8.7) 

2.2 
(8.7-9.9) 

8.02 
(9.9-15.5) 

22.6 
(5.0-15.5) 


5.85 
5.85 
5.85 
6.24 


1.6±0.5[33] 
1.0lt°;£[32] 

1.9±°;I[29] 


7.3 


« 15.6[29] 
(6.7-8.4) 


5(i Ni 
Figs. 3,8 


6.7 


0.24 


10.02 


5.5 
(8.5-11.0) 
5.94 
(11-15) 
11.68 
(6-15) 




78 7Vi 
Figs. 4,9 


6.12 


0.64 


9.3 


9.5 
(8.5-10.2) 

7.6 
(10.2-15.0) 
17.74 
(5.5-15) 



Figure Captions 



Fig. 1. The Ml strength function (in MeV) for 48 Ca with 

an averaging parameter A = 100 keV. 
Fig. 2. Same as in Fig. 1 but for 208 Pb. 
Fig. 3. Same as in Fig. 1 but for 56 Ni. 
Fig. 4. Same as in Fig. 1 but for 78 Ni. 
Fig. 5. Same as in Fig. 1 but for 100 Sn. 
Fig. 6. Same as in Fig. 1 but for 132 Sn. 
Fig. 7. The Ml strength function for 48 Ca, A = 10 keV. 
Fig. 8. The Ml strength function for 56 Ni, A = 20 keV. 
Fig. 9. The Ml strength function for 78 Ni, A = 20 keV. 
Fig. 10. The Ml strength function for 132 Sn, A = 250 keV. 



